Phase diagram of UGe 2 . 

Whether there are quantum phase transitions ? 
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The phase diagram of several itinerant ferromagnets reveals the common feature. The phase 
transition temperature decreases with pressure increase and reaches zero value at some critical 
pressure P c such that at low enough temperatures one can expect critical behavior specific for 
quantum phase transition. It is not the case, however. Being the second order at ambient pressure 
the transition from paramagnetic to ferromagnetic state at high pressures - low temperatures is 
transformed to the discontinuous jump. 

We discuss the magneto-elastic mechanism of development of the first order type instability at 
the phase transition to the ferromagnet state in strongly anisotropic ferromagnet UGe2- Using 
the parameters characterizing the properties of UGe2 we argue the effectiveness of this mechanism 
transforming the very weak first order type transition to the really observable one. 
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I. INTRODUCTION 



The pressure-temperature phase diagrams of several weak ferromagnets exhibit similarity. The transition from the 
paramagnetic to the ferromagnetic states at ambient pressure occurs by means of the second order phase transition. 
The phase transition temperature decreases with pressure increase such that it reaches zero value at some critical 
pressure P c . At some pressure interval below P c the ordered ferromagnetic moment disappears discontinuously. Thus 
at high pressures and low temperatures the ferromagnetic and paramagnetic states are divided by the first order type 
transition whereas at higher temperatures and lower pressures this transition is of the second order. Such type of 
behavior is typical for MnSi [1-4], itinerant ferromagnet-superconductor UGe2 [5,6] (see Fig.l), ZrZn 2 [7]. The same 
behavior has been established in the ferromagnetic compounds Co(Sii_a;Se x )2 [8] and (Sri_ 2; Ca a ;)Ru03 [4] where the 
role of governing parameter plays the concentration of Se and Ca correspondingly. Also, there was demonstrated clear 
evidence for the first order nature of the ferromagnetic transitions in typical ferromagnets like Ni, Fe and Co [9]. 
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FIG. 1. The schematic phase 
diagram of UGe2 Ref. [10]. Thick lines 
denote first order transitions and fine lines 
second order transitions. The dashed line 
is a crossover. Dots mark the positions of 
critical (tricritical) points. 



Here we discuss the magneto-elastic mechanism of development of the first order type instability. Actually the mean 
field treatment of the magneto-elastic mechanism has been put forward in the paper [11] where it was demonstrated 
that the change of transition character from the second to the first order takes place at strong enough steepness of 
the exchange interaction dependence on interatomic distance and large compressibility. It can be considered in frame 
of the Landau theory of the phase transition. Namely in neglect the shear deformation the free energy density near 
the phase transition to the Ising type ferromagnet has the following form 



F = ao(T — T C )M 2 + (3M 4 



K 



qeM 2 . 



(1) 



Here, M is the magnetization density, e is the relative volume change, K is the bulk modulus. The coefficient q is 
related to the Curie temperature pressure dependence as 



a - 



-a K 



de ~ u " dP ' 

In stress absence §f = 0, the deformation is determined by square of magnetization e = j^M 2 that yields 



F = aM 2 + (3 



2K 



M . 



(2) 



(3) 



Hence, at ^ > f3 the phase transition changes its character from the second to the first order. This inequality can 
be rewritten through the measurable parameters as 



KAC 



dT. 
dP 



> 1, 



(4) 



where we used the formula AC 



2^ T c for the specific heat jump at phase transition of the second order. 
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The magneto-elastic interaction also produces another general mechanism for instability of second order phase 
transition toward to the discontinuous formation of ferromagnetic state from the paramagnetic one. For the first 
time it was pointed out by O. K. Rice [12] who has demonstrated that at small enough distance from the volume 
dependent critical temperature T C {V), where the specific heat C/;(r) <~ t~ q , t = T J V) — 1, tends to infinity due to 

the critical fluctuations, the system bulk modulus K = — V 
F = F + Ffi, Ffi <~ -T c T 2 ~ a starts to be negative 



dP 

dv 



V 



d A FV 
dV 2 ' 



T C (V) 

expressed through the free energy density 



K = K - A 



C fl (r)V 2 fm 
dV 



= Kn- AK, 



2 C f l(T) {die 

dP 



< 



(5) 



T->0 



that contradicts to thermodynamic stability of the system. In reality, before there will be reached the temperature 
corresponding to K = the system undergoes the first order transition, such that to jump over the instability region 
directly in the ferromagnetic state with finite magnetization and related to it striction deformation. This transition 
is similar to the jump over the region with dP/dV > on the van der Waals isotherm at the liquid-gas transition. 
The condition of the first order instability (5) can be written in similar to Eqn.(4) form 



K C fl (T) /VTA' 

t c \dp) 



(6) 



However, unlike to Eq. (4) this formula demonstrates that the first order instability is inevitable due to infinite 
increase of fluctuation specific heat. Thus, if in the system with the fixed volume the phase transition is of the 
second order with the infinite increase of specific heat then the effect of finite compressibility under assumption 
that the critical temperature is the volume dependent parameter transforms it into the phase transition of the first 
order. In reality, the striction interaction can change the shape of the free energy singularity in respect to its form 
in incompressible case. More elaborate treatment [13] taking into account this effect leads to the following condition 

of the first order instability jr 3 ^^ M f"( x ) (tjp ) > 1- Here the function f(x) determines the fluctuation part of free 



energy 



F 



-T, 



T-T c 



fi is the shear modulus. Usually, the left hand side in Eqn. (6) is quite small and the 

transition of the first order occurs at temperature T* close to the critical temperature where fluctuation specific heat 
is large enough. It means that the temperature difference T* — T c is smaller than the critical temperature T c by many 
orders. The latent heat at this transition q w Cfi(T* — T c ) proves to be extremely small. So, the first order phase 
transition is practically indistinguishable from the second order one and called weak first order phase transition or 
the phase transition of the first order closed to the second order. 

According to Eqs (4), (6) the magneto-elastic mechanism effectively leads to the first order transition when the 
critical temperature is strongly pressure dependent. This is the case in all mentioned above materials. To check 
the criteria (4), (6) one must calculate the mean field jump and fluctuation part of the specific heat near the Curie 
temperature for each given material. To be concrete, here, I'll do these calculations for UGe2 characterized by strong 
magnetic anisotropy and by the precipitous drop of the critical temperature at pressure increase near 14-15 kbar [14]. 



II. THE SPECIFIC HEAT NEAR THE CURIE TEMPERATURE 



UGe2 is orthorhombic crystal with ferromagnetic order at ambient pressure found below T c = 53 K. Magnetic 
measurements reveal a very strong magnetocrystallinc anisotropy [15] with a being the easy axis. We shall denote it 
as z direction. The free energy of strongly anisotropic ferromagnet can be written in terms of one component scalar 
order parameter corresponding to magnetization density M z (r) along z axis. In that follows we shall omit the order 
parameter index z. 



(7) 



Here, a — a (T — T c ), the gradient terms are written taking into account the orthorhombic anisotropy 7^ = 
7xx \ 

j yy , where the x, y, z are directions of the spin axes pinned to b, c, a crystallographic directions corre- 
lyy ) 

spondingly. The last nonlocal term in Eq. (7) corresponds to magnetostatic energy [16,17] — MH — H 2 /8ir, where 
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internal magnetic field H expressed in terms of magnetization density by means of Maxwell equations rotR = and 
div(H + 47rM) = 0. We shall use the following estimations for the coefficients in the Landau free energy functional 

a a = — (8) 
m z n 

P= , T 2 \ 2 , (9) 

7-«7y«7.«^- (10) 

Here, m = 1.4/is is the magnetic moment per uranium atom at zero temperature [18], n — a -3 is the density 
of uranium atoms, which can be approximately taken equal to inverse cube of the nearest-neighbor uranium atoms 
separation a — 3.85 Angstrom [19]. 

The mean field magnetization and the jump of specific heat are 

= »■ < 12 > 

The experimentally found specific heat jump AC exp w 10 m ^ lK ~ 1 per uranium atom [19] is in remarkable corre- 
spondence with Eq.(12). 

The effective Hamiltonian of noninteracting field of the order parameter fluctuations is given by 

H = J2( a + Hjhkj + 2nk 2 /k 2 ) M k M_ k , (13) 

k 

where Mk = / M(r)e~ lkr d 3 r. The corresponding free energy and the specific heat are [20] 



_ T 2 al r dk x dk y dk z 
fL0 ~ 2(2^ J [a + 2Kkl+ ll3 k t ktf 1 ' 



Proceeding to spherical coordinates and performing integration over modulus k we come to 



r T 2 a 2 r 1 



<K l (a + 27rC 2 )V2( 7± + C( lz - 7± ))3/2 • ( 16 ) 



Here, ■~f±((p) — j x cos 2 tp + j y sin 2 <p. At critical temperature a = and the integral diverges. Hence, performing 
integration over £ with logarithmic accuracy we obtain 



Cno = In — w \/ ^-ln , (17) 

} 32tt\/2^ : 7 3 / 2 2tt 32tt V 27rm 2 n T — T c v ; 



where 



7 3/2 2tt7 7 3/2 (yj) - 

The used condition a -C 27r at T c = 10AT is realized at 

T-T c 27rm 2 n 

< « 0.015. (18) 

c J c 

In view of roughness of the parameter estimation the region of logarithmic increase of specific heat can be in fact 
broader. 
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The calculation taking into account the interaction of fluctuations has been performed by Larkin and Khmelnitskii 
[21]. In our notations the expression for the fluctuation specific heat at const pressure obtained in this paper is 

Ql/3 J-2 a 2 1/3 

C f i= 'f" (lnf ) (19) 

Here "/lk = ^/Jfi'V? 1S ^ nc effective constant of interaction. Using the Eqs. (8)-(10) one can rewrite Eq. (191 as 

V6 / o^2„ \ 1/3 



^"™U==kr ■ (20) 



10 \2irm 2 n J \ T — T c 

So, we see that the order parameter fluctuations give rise the increase of specific heat near the critical point. The 
power of the logarithm (In ^r) 1 ^ 3 is quite slow function slightly exceeding unity, hence in the temperature region given 
by inequality (18) one may estimate the fluctuation specific heat as 

n 

Cfi > 5. (21) 

We see that the fluctuation specific heat is smaller than the mean field jump given by Eqn. (12). Hence to check the 
first order phase transition instability in UGe2 one must to proceed with criterium (4). 

III. INSTABILITY OF THE SECOND ORDER PHASE TRANSITION 

The Curie temperature in UGe2 falls monotonically with increasing pressure from 53 K at ambient pressure and 
drops precipitously above 15 Kbar [14]. The average value of the critical temperature derivative can be estimated as 

dT c 40 Kelvin „ ^ q 

— £ « = 4 x 10~ 2j cm 3 (22) 

dP 14 kbar v ; 

For the bulk modulus we have 

K = pc 2 « 10 n erg/cm 3 , (23) 

where we have substituted typical sound velocity c ss 10 5 cm/ sec and used known [22] density value p — 10.26 g/cm 3 . 
Thus, we have for the combination Eq.(4) 



Kn (dT c 



= 0.2 . (24) 



At T 1=3 \QK the pressure derivative of the critical temperature is much higher (and its square is even more higher) 
than its average value given by Eq. (22). So, we come to conclusion that at critical temperature of the order 10 K 
the criterium (4) is fulfilled and the phase transition of the second order turns into the first order one. 

IV. CONCLUSION 

The magneto-elastic interaction provides development of the first order instability at the phase transition to the 
ordered state in a ferromagnet. However, actual temperature interval of this instability development is negligibly 
small and the first order transition looks almost indistinguishable from the second order one. The particular feature 
of anisotropic ferromagnet UGe 2 is the precipitous drop of the Curie temperature as the function of pressure near 
14-15 kbar. Due to this property at about these pressures the second order phase transition (or very weak transition 
of the first order) to ferromagnet state turns into the real first order type transition. 

At low temperatures according to the Nernst law and the Clausius-Clapeyron relation 



dT c vi - v 2 



dP S\ — S2 



-> 00 (25) 



the drop of transition temperature with pressure begins to be infinitely fast. It means that weak first order transition 
has the tendency to be stronger and stronger as temperature decreases. Hence, the effect of magneto-elastic interaction 
or, more generally, of the order parameter interaction with elastic degrees of freedom at arbitrary type of ordering 
raises the doubts upon the existence of quantum critical phenomena. 
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